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APPROXIMATE AERODYNAMIC INFLUENCE COEFFICIENTS FOR 
WINGS OF ARBITRARY PLAN FORM IN SUBSONIC FLOW 
By Franklin W. Die derich 

SUMMARY 


Aerodynamic influence coefficients for symmetrically loaded wings 
of arbitrary plan form in subsonic flow are derived from a simple 
empirical method of estimating spanwise lift distributions. The 
application of the coefficients to an aeroelastic analysis is 
discussed. 


INTRODUCTION 


In the aeroelastic analysis of reference 1 the spanwise lift dis- 
tribution was assumed to be given by strip theory with over-all reduc- 
tion and rounding at the tip to take aerodynamic induction approxi- 
mately into account. The spanwise lift distribution may be incorporated 
into the aeroelastic analysis more accurately by using suitable aero- 
dynamic influence coefficients. In this paper a method is presented 
for determining approximate aerodynamic influence coefficients for 
symmetric loadings from the method of estimating lift distributions 
presented in reference 2. An approximate means of taking the effect 
of the fuselage Into account is also presented in the present paper; 
it is based on the assumption that the fuselage carries an amount of 
lift which depends solely on the angle of attack of the wing at its 
root. 


Although this paper is intended primarily as an adjunct to refer- 
ence 1, the method and coefficients presented herein may be used for 
calculating lift distributions apart from reference 1. In order to 
facilitate the application of the present paper to such calculations, 
all required definitions and numerical factors are presented herein. 

As in reference 1 matrix notation is used in the present paper. A 
short summary of matrix notation and algebra and a discussion and 
derivation of the pertinent integrating matrices are given in 
reference 1. 
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SYMBOLS 


b 

b* 


c 

c 

c e 

c s 

c l 

C La 

w 

f 

g 

L 1 ^ 

H 

s 

y 

y* 

a 

a 

T1 

T]* 


span 

span less fuselage width 

chord, parallel to plane of symmetry 

average chord ^ j 

average of root chord and chord in plane of symmetry 

chord in plane of symmetry 

local lift coefficient 

wing lift -curve slope 

fuselage width 

"ideal" function (see fig. 1) 
factor defined by equation (11) 

integrating matrices 
first row of Kq_ matrix 


aerodynamic-influence-coefficient matrix 
total wing area 

lateral ordinate measured from plane of symmetry 

dimensionless lateral ordinate 

angle of attack 
average angle of attack, defined by equation ( 2 ) 
lateral ordinate measured from wing root 

dimensionless lateral ordinate 
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^ taper ratio 

A angle of sweepback at quarter -chord line 

Subscripts: 

a additional 

t basic 

r root 

Matrix notation: 

square matrix 
diagonal matrix 
column matrix 
row matrix 
■unit matrix 

CALCULATION OF INFLUENCE COEFFICIENTS 
Wing Alone 

In reference 2 the spanwise lift distribution at subsonic speeds 
for any angle -of -attack distribution is obtained by considering sepa- 
rately the additional lift distribution and the basic lift distribution. 
The additional lift distribution is calculated from the average of the 
chord distribution and an "ideal" function f (given in reference 2 and 
reproduced in fig. 1 of the present paper), which has been obtained 
empirically for various angles of sweep: 



(tt)* " i(t + 


(1) 
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The average angle of attack a is obtained from 

dy* 


a 


2 r b / 2 

* B J 0 “ C * 


= T 1 a-£ 
j 0 c 


(2) 


where y* is the dimensionless lateral ordinate defined by 


y* = 


hj2 


Similarly, the basic lift distribution is obtained from one -half 
of the product of the chord distribution and the angle -of -attack 
distribution for zero lift (a - a) : 

(^- **!*“-*> (3) 

Any sharp corners in the resulting distribution must be rounded off. 

Equations (l) and ( 3 ) may be combined and written in matrix 
notation as 





w 


where |lj represents a column of l*s. The integral required in cal- 
culating a may be evaluated by means of an Integrating matrix. The 
first row of the Kq matrix' of reference 1 IjCqJ serves to perform 
this operation and round off the basic distribution at the same time . 
The six-point jjCqJ matrix (shown in table 1 of the present 

paper) is based on spanwise stations located at 0, 0.2, 0.4, 0.6, 0.8, 
and 0.9 of the distance over which the integration is performed. Con- 
sequently, the analysis of the wing alone is based on stations 
corresponding to y* = 0, 0.2, 0.4, 0.6, 0.8, and 0.9» 

In order to fulfill the requirement that the integral of the basic 
lift distribution is zero, the value of a may be determined from the 
relation 
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where, for the special case of the wing alone. 



(6) 


The value of g is presented for various values of the taper ratio in 

figure 2. The case of the wing alone corresponds to — = 0. 

"b 


Equation (5) may he substituted in equation (4) and both sides 

multiplied by — to yield 
c r 



( 7 ) 


where the aerodynamic-influence-coefficient matrix Q is defined by 



(8) 


The Kq matrix is a square matrix consisting of rows all equal to jjCqj^ 
The elements of ^|lcq and. ^[jeJ are given in table 1, the factor 


of being included to .avoid the necessity of a final multiplication by 

that factor. In equations (7) and (8) the chord at the plane of sym- 
metry could be used instead of the root chord c r since the two are 
identical for a wing without fuselage. 


In calculating the Q matrix the f values are first read for 
the given sweep angle from figure 1 at y* values of 0, 0.2, 0.4, 0.6, 
0.8, and 0.9} these f values are then entered in a diagonal matrix as 
shown in table 2 for an example wing with 30° sweepback. A value of g 
is then read from figure 2 for the given taper ratio (X = 0.5 for the 
example wing) and multiplied into the f matrix. -The resulting matrix 

is postmultiplied by the ||jtqJ matrix of table 1. (See table 2 for 
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the result of this operation in the case of the example wing. ) The 
matrix is then added to the product. The resulting matrix is 

Q 

postmult iplied by a diagonal matrix of the — values at the previ- 

c r 

ously cited values of y* in order to obtain the desired Q matrix. 
(See table 2.) 


Wing and Fuselage 

For the purpose of calculating aeroelastic effects the presence of 
the fuselage may be taken into account in several ways. In one method 
the wing root is considered to be at the plane of symmetry and the part of 
the wing inside the fuselage, to be infinitely stiff. In this case the 
coefficients obtained in the foregoing section may be used. 

Another method consists of calculating the aerodynamic influence 
coefficients for stations corresponding to tj* (rather than y*) values 
of 0, 0.2, 0.4, 0.6, 0.8, and 0.9- In that case an assumption has to be 
made concerning the amount of lift carried by the fuselage and the 
effect of this lift on that carried by the wing. 

For untwisted wings it is commonly assumed that the fuselage car- 
ries the same lift as the part of the wing it intercepts. In this 
analysis the same assumption is made for a twisted wing, and the part of 
the wing intercepted by the fuselage is assumed to be at the angle of 
attack of the wing root. The value of a for this assumption may then 
be obtained from. the requirement that the basic lift distribution 
integrates to zero: 



where w, b’, c s , and c r are geometric parameters defined in 
figure 3, 


Cp = 


c s + c r 


( 10 ) 

; 

w 


O 
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g = 


£t+NJ]§ 


( 11 ) 


°5 + S 8Xe P resented - In figure 2 as a function of the fuselage- 
width and' taper ratios. 6 

.. _y ±th the value °1 a given hy equation (9). equation (k) takes 
the form of equation (7) where Q is now defined hy 


H- + l&il [1] 


( 12 ) 


Comparison of equations (9) and (12) with equations (5) and (8), respec- 
tively^ indicates that the Hcpl matrix is obtained by adding the 
w c e L ■ 

quantity — — to every element of the first column of the [f-jl 
matrix. L J 


When the Q matrix of equation ( 12 ) is used in equation (7), the 

values of will be obtained at stations corresponding to tj* = 0 , 

0.2 O.k, 0.6 0.8, and 0.9 for the given values of the angle of attack 
at those stations. The relation acJ£ 


_ y. U w V 


may facilitate determination of the stations of interest. In calcu- 
lating the Q matrix for the wing-fuselage combination, values of the 

f function and the ^ ratio are obtained at those stations. The 

average chord c e , which is calculated from equation ( 9 ), is used to obtain 

the term ~ — which is added to each element of the first column of 


the — 


Its. 
2 


matrix of table 1 to yield the — 

2 


Mfle 


matrix. 


ma ^ iX P rem ultiplied by a diagonal matrix obtained by 
multiplying the f matrix by the value of g taken from figure 2 for 
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the given configuration. The matrix is added to the resulting 

c 

matrix and the result is postmultiplied by the — . matrix to yield 
the desired Q matrix. 

APPLICATION OF INFLUENCE COEFFICIENTS 


The aerodynamic influence coefficients calculated as outlined in 
the foregoing section may be used to calculate the lift distribution of 
wings with arbitrary sweep and taper for any. symmetrical angle -of -attack 
distribution. (For discontinuous distributions, such as those caused by 
flap deflections, the equivalent 8 values of reference 3 may be used 
in order to obtain the rounding of the lift distribution at the point of 
discontinuity called for in reference 2.) Once the coefficients have 
been calculated for the plan form of interest (with or without fuselage), 

CC7 

the values of the loading coefficient at stations corresponding 

i c r 

to 0, 0.2, 0.1(-, 0.6/ 0.8, and 0.9 of the distance from the wing root to 

the wing tip are obtained by premultiplying a column matrix of the 
values of the angle of attack at those stations by the influence - 
coefficient matrix. 


When the influence coefficients are used in conjunction with the 
method of aeroelastic analysis presented in reference 1, steps © and 


in table VI of reference 1 must be modified. 

c 


The influence - 

matrix used in step (l^ . A matrix 
must be calculated ^from the given values of the moment 


coefficient matrix replaces the 
_ o 
e l _c_ 

_ el r Cr _Jsub 


arm 


c \ 

ej_ pertinent to subsonic speeds and the given chord ratio — ) 

Cr / 


and postmultiplied by the influence -coefficient matrix. 

e ! ' ~ ' 2 


matrix replaces the 


l5r 


(t) 


The resulting 

matrix called for in step © . This 


—•sub 

procedure is not applicable to the supersonic case since the aerodynamic 
influence coefficients have been obtained only for "the subsonic case. 
Elsewhere in the method of reference 1, the lift -curve slope CLq. of 

the present paper is used in place of the effective slopes mg and 
mej. 




the ratio 


mg 


is taken as 1, and the equivalent angle of attack a 


of reference 1 (not to be confused with the average angle of attack a 
of the present paper) is merely the sum of the geometric and structural 
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angles of attack. A convenient way of obtaining the lift coefficient 
for any calculated loading consists of multiplying the a values 
obtained from equation (5) or (9) by the lift-curve slope Cl^. In all 
other respects the calculations indicated in reference 1 are unaffected 
by the use of the influence coefficients. 

Compressibility corrections may be applied, as suggested in refer- 
ence 2 , by using an f function corresponding to an effective sweep 
angle and by using a corrected lift-curve slope. 


DISCUSS 100 


The influence -coefficient matrices of this paper constitute an 
expression of the empirical method of reference 2 in matrix form. The 
matrix elements are influence coefficients in the generalized sense that 
if multiplied by the angle -of-attack values along the span and summed up 
they will yield the value of the lift at a given station. Individually 
the elements do not necessarily represent the value of the lift at a given 
station due to a unit angle of attack at another station in the more 
narrow meaning of the term "influence coefficient.” 

The discussion in reference 2 of the merits and demerits of the 
empirical method applies to the aerodynamic influence coefficients as 
well. In short, at subcritical speeds they may be expected to yield 
results which are more reliable, than those of a strip-theory analysis 
but not so reliable on the average as the results of one of the more 
refined theories for spanwise load distribution. The results differ 
from those of the method of reference 2 in only one respect. The 
rounding of the basic loading curve required in reference 2 is performed 
by integrating matrices in the present paper. In order to show the 
effect of this difference, the results of the two approaches for the 
basic lift distribution of an untapered unswept wing with washout of 
1 radian are shown in figure k. Also shown are results calculated by 
means of Multhopp's and Weis singer’s analytical methods and on the basis 
of the modified strip theory used in reference 1 . 

The influence-coefficient method and the method of reference 2 
yield almost identical results, which are in fair agreement with the 
results of Multhopp’s and Weissinger’s analytical methods. The agreement 
of the results calculated by the modified strip theory of reference 1 with 
the analytical results is somewhat poorer on the average. (The results of 
unmodified strip theory would be much poorer; for the case of fig. k the 
loading-coefficient curve would be 80 percent larger at every point t han 
that of the modified strip theory and would not be rounded off at the tip.) 
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The methods used for the calculations represented in figure k 
(except that of modified strip theory) would yield almost equal 
results for the additional loading, which would he intermediate between 
elliptical and rectangular^ whereas the result of modified strip theory 
would be a rectangular distribution rounded only at the tip. The 
results of the influence -coefficient method are therefore superior to 
those of the modified strip theory for both basic and additional lift 
distributions, although the difference is not so large for tapered wings 
as for the untapered-wing considered in figure 4. Since the agreement 
of the results of the influence -coefficient method with those of the 
analytical methods is somewhat poorer for basic lift distributions than 
for additional lift distributions, the influence coefficients may 
be expected to be most useful when the angle of twist is less than the 
angle of attack. The influence coefficients will probably furnish 
fairly reliable results for the change in loading due to aeroelastic 
effects but somewhat less reliable results for the divergence speed. 

In general, they tend to overestimate the lift, particularly at the wing 
tip, so that aeroelastic effects are overestimated to some extent. For 
similar reasons the influence coefficients of this paper- are not appli- 
cable to antisymmetric lift distributions because these lift distri- 
butions do not have an additional part. 


CONCLUDING REMARKS 


A method has been presented for calculating aerodynamic influence 
coefficients for symmetrical loadings at subsonic speeds from a simple 
empirical method of estimating spanwise lift distributions. The coeffi- 
cients axe particularly suited for aeroelastic analyses and may be 
expected to give good results for the changes in spanwise loading 
associated with aeroelastic effects but probably less reliable results 
for the divergence .speed. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., March 3, 1950 
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Figure 1.- The "ideal" distribution f 
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Figure 2.- Values of the factor g 
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Method 

— Influence coefficients 

— Reference I 
--Referenced 

— Mu Ithopp 

— Wetssinger 


Di mensionless lateral ordinate, y * . 


Figure 4.- Basic lift distribution of rectangular wing of aspect 
ratio 5 with unit linear symmetric twist in incompressible 
flow. 
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